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A mean field calculation of the Hubbard model finds a 
rich phase diagram. The antiferromagnetic phase is generally 
unstable away from half filling, and there are several regions of 
phase separation. One solution in particular closely resembles 
the stripe phase of White and Scalapino. By comparison to 
unrestricted Hartree-Fock calculations (for which this phase 
is metastable), it is demonstrated that this phase arises from 
phase separation. The interface surface tension is found to 
change sign below a particular stripe width, at which point the 
stripes begin to meander, gradually crossing over to diagonal. 

When the Hubbard model is doped away from the an- 
tiferromagnetic insulator at half filling, a number of cal- 
culations find evidence for spatially inhomogeneous solu- 
tions. There is considerable debate as to whether 
these solutions are generic features of the Hubbard 
model, or arise only in a restricted parameter domain. 
Related issues are whether the inhomogeneity is driven 
by phase separation or antiferromagnetic (AFM) domain 
wall formation ||^, and how these features are related 
to 'stripes' in cuprates and other oxides. Similarly, in 
the doped tJ model, various calculations find that the 
ground state is striped or uniform or phase sep- 
arated 1^,0. Inclusion of realistic values of t' into the 
model further reduces stripe stability Reference to 
earlier calculations may be found in these articles and in 
the reviews 

Unrestricted Hartree-Fock (UHF) calculations 10-1^] 
find that the holes form filled (one additional hole per 
row) stripes which act as antiphase boundaries between 
AFM domains. Such filled domain wall stripes are not 
found in more advanced calculations [^,|| of the Hub- 
bard model, and are not consistent with experiment on 
the cuprates We here analyze a metastable state 

of the UHF calculations, which closely resembles the 
White-Scalapino (WS) stripes, and agrees better with 
experiment. These stripes can be understood from a 
phase separation approach, comparing the free energies 
of low-order commensurate magnetic phases, qx,<ly ~ 
or Qi — ■n/a. The resulting mean-field phase diagrams 
involve phase separation between the AFM phase and 
a metallic phase, either ferromagnetic (FM), as in early 
ferron phase approaches to the Hubbard model , or a 
phase resembling WS stripes, depending on the value of 
second neighbor hopping parameter t' . These stripes are 
stable local free energy minima in UHF calculations, but 
globally there are alternative states of lower free energy 
[ p^ . However, these solutions can be stabilized by addi- 
tional interactions beyond the pure Hubbard model (e.g., 
charge-density wave or superconducting) , and hence may 
be relevant to experiment. These additional interactions 



will be discussed in a companion publication [|l5| ; here 
we introduce the mean-field model and utilize UHF cal- 
culations to calculate the surface tension in the resulting 
stripe phases. We find that WS-like stripes are stable 
against macroscopic phase separation. 

We study a one-band electron-hole symmetric Hub- 
bard model [interaction = UYliii^i] — l/2)(ni| — 1/2)] 
with bare dispersion tk — —2t{cx + Cy) ~ At'cxCy, with 
Ci = cos kiO. In the presence of a mean-field magneti- 
zation m, 
becomes 



q at wave vector q, the quasiparticle dispersion 
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E± - -(efe-|-efe+, ±£;o) 



(1) 



where Eq = y (^^ ^ ^fc+g)^ + 'iU^m'^. The site magneti- 
zation is found self-consistently from 



(2) 



with Fermi function f{E) = 1/(1 + Q(E-E,.)/kBTy r^-^e 
free energy is 



k,i=± 



(3) 



where S is the entropy. 

The phase with q ~ Q = (tt, tt) provides a good model 
for the AFM phase at half filling with Mott gap, success- 
fully describing the spin wave dispersion |16| and Monte 
Carlo results jljjl^, and serving as the basis for a num- 
ber of treatments of strong correlation effects [|9|. A fit 
to the dispersion of the magnetic insulator SrCu02Cl2 
finds |§ t = 325meT/, U = 6.03t, and t' = -0.276t, 
which will be assumed below unless stated otherwise. We 
use the same model, with different choices ol q ^ Q, 
to describe a number of competing magnetically ordered 
states. While at half filling the AFM state has lowest 
free energy, this is not true for finite doping x, leading 
to a rich phase diagram, with regimes of phase separa- 
tion. Figure [l] shows the low-temperature free energy 
as a function of doping for the case i' = for three 
magnetic phases, the standard antiferromagnet (AFM) 
q^Q,a, ferromagnet (FM) with q — (0, 0), and a linear 
antiferromagnet (LAF) with q — (tt, 0) (see Fig. be- 
low). The curves are symmetric about half filling (x = 0). 
For \x\ < 0.25 the AFM lies lowest in energy; between 
0.25 < \x\ < 0.65 the LAF lies lowest, and beyond that, 
the ground state is nonmagnetic {rriq = 0). For all dop- 
ings, the FM state is metastable. At high doping the 
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magnetic phases terminate when — > 0. The inset to 
Fig. d shows the dispersions for the stable phases, AFM 
at x = (sohd hnes) and LAF a,t x — 0.353 (dashed 
lines). 




FIG. 1. Free energy vs. doping for several magnetic phases 
of the Hubbard model {U = 6.03t, t' — 0). Diamonds = 
AFM, triangles — LAF, circles — FM, and squares — PM 
phase. Dashed lines = tangent construction. Inset: Disper- 
sion of magnetic phases: solid lines — AFM at x = 0, dashed 
lines = LAF at a; = 0.353; Brillouin zone points F = (0,0), 
X = (^,0), S = (7r,^). 

The antiferromagnetic state has a cusp at half fill- 
ing, with the slope discontinuity being the Mott-Hubbard 
gap. Away from half-filling, this state is always thermo- 
dynamically unstable - the compressibility d'^f/dx^ 
is negative. The tangent construction (dashed lines) 
shows that the equilibrium state between zero doping 
and |a;| — Xc = 0.353 consists of a phase separation be- 
tween the AFM and LAF phases. Note that the mean- 
field model misses the true UHF ground state, which has 
filled {x = 1) stripes in an AFM background. It can 
be shown that, if the last term in Eq. ^ is omitted, this 
ground state would be recovered for large U, with the 
LAF phase stable only for a small parameter range near 
U/t = 6-8 (f = 0). 

The resulting phase diagram x vs U is shown in Fig. |^. 
Phase separation persists for all finite U, but while the 
insulating state is always AFM, there is a crossover in 
the metallic stripe component from paramagnetic phase 
for U <Uc = 5.3t to LAF for U > Uc). When t' ^ 0, the 
phase diagram is completely different, with phase sepa- 
ration between the AFM and a FM phase . For large 
U, the Hubbard model should reduce to the tJ model; 
agreement with recent calculations for the phase separa- 
tion boundary in the Hubbard [|j and tJ Q models is 
satisfactory (triangles and -l-'s in Fig. ^a). The deviation 
at small U (large J) is expected, since the models are 
equivalent only in the large-C/ limit. While the metallic 
phase in the tJ model is usually taken as paramagnetic, 
the WS results may hint that it is an LAF phase near 
U = lit. 
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FIG. 2. Phase diagram, x{U) for the Hubbard model, with 
t' = 0. Triangles = Hubbard model results, estimated from 
Fig. 1 of Ref. PI; dot-dashed line with -f 's = tJ model results, 
Ref. |]b, assuming J/t = At/U. 

Figure ^ illustrates some of the low-energy textures 
found in UHF calculations, and shows that in the LAF 
stripe phase dispersions the added states form addi- 
tional bands near midgap, as found in ordered stripe 
arrays pO[ | and for randomly distributed magnetic po- 
larons (Fig. ||a,b). [For the left-hand panels of Fig. ||, 
the UHF calculations were iterated to self-consistency on 
24x24 (a), 32x6 (c), or 12x12 (e,g) lattices with periodic 
boundary conditions. For the dispersions of the right- 
hand panels of Fig. ^ these solutions were extended to 
a 32x32 (d,h) or 36x36 (f) lattice, with one additional 
iteration (Fig. ||b was on a 24x24 lattice).] 

The LAF stripes resemble the WS stripes of the tJ 
model: the minimum LAF stripe is two cells wide, and 
acts as an antiphase boundary between AFM domains. 
Fig. ||c. In both calculations, the doped ground state 
is found to involve mixtures of LAF and AFM stripes, 
with no sign of insulating, empty stripes. The doping is 
comparable: the star in Fig. ||a represents the WS stripe, 
assuming an effective U/t = 4:t/J, with J = 0.35i. Both 
kinds of stripe have similar fractional transfer of holes 
onto adjacent AFM rows (see caption of Fig. and both 
are destabilized by non-zero t'. We find that the charged 
stripes have a fixed, minimal width for x < 1/6, with the 
charge per row of a stripe doubling at higher doping, and 
the stripe phase terminating near x — 1/3; WS find simi- 
lar doping dependences, systematically shifted due to the 
difference in hole density (1/3 vs 1/4) on a stripe. Simi- 
lar LAF stripes were found previously as metastable UHF 



solutions |11|. An LAF- like state has also been found in 
recent Monte Carlo calculations in the manganites [ pl| ; 
interestingly, a spin fiux phase can form from a coherent 
superposition of two LAF phases. 
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FIG. 3. Unrestricted Hartree-Fock calculations of the Hub- 
bard model. Left panels: the hole configurations for 1/8, 1/6, 
2/9, and 1/3 dopings minimizing the free energy. Right pan- 
els: the corresponding dispersions (open circles); solid lines = 
mean-field bands for undoped antiferromagnet (b) and LAF 
(h); dashed lines = chemical potential. (a,b): Random distri- 
bution of holes in an antiferromagnetic background forming 
magnetic polarons at 1/8 doping. (c,d): LAF stripes at 1/6 
doping - on the charge stipes the hole density is 0.27 and the 
magnetization is 0.46; on the antiferromagnetic stripes the 
hole density is 0.06 and the magnetization is 0.77. (e,f): AFM 
stripes in LAF background at x=2/9; there is a weak mod- 
ulation of the hole density on the LAF stripe: 0.245 (0.365) 
holes on the outer (inner) rows. (g,h): LAF configuration at 
1/3 doping. 

It is clear from Fig. || that the LAF- AFM stripes arise 
from phase separation, with the stripe spacing evolving 
as expected with doping. However the question remains 
as to whether the stripes are stable against macroscopic 
phase separation - i.e., is the free energy of the stripe 
phase Fstr higher or lower than that of the separated 
bulk phases (tie- line) i^sep? We answer this via UHF cal- 
culations, taking care to minimize finite size effects. This 
is done by (a) adjusting U so that Xq is a simple ratio- 
nal fraction (at U/t = 8, — 1/3), and (b) working 
with large lattices, up to 128 x 6. Figure ^ shows (a) 
the average free energy on each row of a series of AFM- 
LAF stripe arrays, of the same average doping [x = 1/6) 
but different stripe widths and (b) the resulting surface 
tension a (free energy difference Fstr — F^ep per domain 
wall atom). The magnetic contribution, associated with 
excess holes pushed onto the magnetic bounding layers, 



is positive and saturates for wider stripes, while the LAF 
energy is negative, and oscillates on the stripes, in par- 
allel with hole density oscillations. These oscillations are 
due to quantum confinement, similar to the Friedel os- 
cillations seen in electrons confined on a step on a Cu 
surface [|2^ . The confinement oscillations lead to a long- 
range interaction between domain walls (across a charged 
stripe) , which explains why the net surface tension satu- 
rates so slowly as a function of stripe width, and why it 
depends mainly on the width of the LAF stripes. 
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FIG. 4. (a) Free energy f{x) per row for a LAF stripe array 
in the Hubbard model, with t' = 0, U = St and a;=l/6, 
comparing several different lattice periodicities: 32 LAF x 32 
AFM for a 128x6 lattice, and, for 64x6 lattices, 16x16, 8x8, 
and 2x2. Energies of narrower arrays are offset for clarity, (b) 
Excess free energy per domain wall atom Edw plotted vs 1/A'', 
where A'' is the charge stripe periodicity, for a; = 1/12 (circles), 
1/9 (squares), 1/6 (diamonds), and 2/9 (stars). Dotted lines 
connect unstable vertical stripes, (c) Comparison of energy 
per site Esuc of meandering (circles) and diagonal (squares) 
stripes compared to the filled diagonal stripes (diamonds) . 
(d) Pattern of spin and charge order on meandering a;=l/6 
stripe. 

For the widest stripes, the surface tension a starts to 
level off to a value of Q.lt per domain wall atom for an 
isolated domain wall. As the stripes move closer a de- 
creases, ultimately changing sign (negative surface ten- 
sion). When the LAF stripe has a width of 8 cells, the 
surface tension is essentially zero. For narrower LAF 
stripes straight vertical stripes are unstable, but can be 
pinned by commensurability effects on specially chosen 
lattices; the free energy is generally high (dotted lines 
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in Fig. IJb). On larger lattices, the UHF spontaneously 
evolves to a meandering stripe pattern, Fig. with free 
energy lower than the tieline (Fig. lowest points of 
solid lines). The meandering LAF stripes are composed 
of straight diagonal segments, separated by kinks. On 
the straight segments, holes on successive rows are shifted 
diagonally by one Cu site, leading to ferromagnetic align- 
ment between kinks. Remarkably, the free energy of the 
meandering stripes is lower than that of the correspond- 
ing straight diagonal stripes. Fig. ^ (although the differ- 
ence is very small, and could be a finite size effect). The 
crossover appears to be kinetic energy driven: the holes 
in the LAF phase are delocalized along the (FM) rows, 
but when the LAF stripes get too narrow, adjacent rows 
shift to provide a FM coupling. Due to commensurability 
pinning effects, it will be hard to repeat this calculation 
for arbitrary values of U , although xq ~ 1/4 at U — l&t. 

There is a gradual crossover (Fig. ||c) from vertical 
(at x=2/9) to meandering (1/6, 1/9) to diagonal stripes 
(1/12). At x=l/12, the diagonal stripes have a low free 
energy, and meandering configurations are unstable. Fig- 
ure l^c also includes the free energy of the diagonal, one 
hole per row stripes which are the UHF ground state 
[ p2t . The free energy differences are small, and the order 
of states may be reversed by including some additional 
(perhaps phononic or Coulomb) interactions. 

Some mention must be made about the size of the lat- 
tice used. Most of the results correspond to lattices 96 x 6 
(for X = 1/9 and 2/9), 128 x 6 (for all x = 1/12, and for 
the largest period at x = 1/6), or 64x6 (for the remaining 
X — 1/6), with periodic boundary conditions assumed. 
The meandering stripes were all on 48 x 12 lattices, and 
the diagonal on 24 X 24 (48x16 for 1/12). Straight stripes 
are metastable when the LAF stripe width is 4, and we 
had to use special lattices to stabilize this configuration: 
12 X 24 for X = 2/9, and N' x 12, with N' = 48 (1/6), 24 
(1/9), and 32 (1/12) (there was no similar problem for 
the LAF width=2 stripes, which are also metastable). 
The surface tension also depends sensitively on the free 
energies of the reference end phases. These can be cal- 
culated either exactly, from the mean field theory, or nu- 
merically from the UHF. For the LAF the agreement is 
quite good: energy per site ELAp/t — 1.46579 (mean 
field) vs 1.46578 (UHF); EAFivi/t = 2.46577 (mean field) 
vs 2.46588 (UHF) (UHF's on 24 x 24 matrices). It was 
necessary to use the UHF value for Eafm to calculate 
surface tensions. 

The present calculations shed some light on the con- 
troversy in the tJ model. The doped AFM phase is so 
unstable in the Hubbard model, that it is likely that 
the elementary excitations in the 'uniform' lightly-doped 
tJ model are really magnetic polarons. A recent Quan- 
tum Monte Carlo study of the Hubbard model ||l^ also 
finds that holes add new dispersionless bands, and do 
not uniformly dope the AFM phase. (See also Ref. [p3[.) 
Hence, the three-sided debate about 'uniform' (or mag- 



netic polaron) ||] vs stripe Q vs (macroscopically) phase- 
separated tJ ground state is in all probability really 
a debate about three kinds of phase separated ground 
state. Our results favor (meandering) stripes. 

In conclusion, we find WS-like stripes at the HF level 
in the Hubbard model (albeit as metastable states), and 
we demonstrate that they arise from a tendency to phase 
separation, providing the first estimate of their surface 
tension. 
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